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ABSTRACT 

The tidal disruption of a star by a supermassive black hole (SMBH) is a highly 
energetic event with consequences dependent on the degree to which the star plunges 
inside the SMBH's tidal sphere. We introduce a new analytic model for tidal disrup- 
tion events (TDEs) to analyze the dependence of these events on /3, the ratio of the 
tidal radius to the orbital pericenter. We find, contrary to most previous work, that 
the spread in debris energy for a TDE is largely constant for all /3. This result has 
important consequences for optical transient searches targeting TDEs, which we dis- 
cuss. We quantify leading-order general relativistic corrections to this spread in energy 
and find that they are small. We also examine the role of stellar spin, and find that 
a combination of spin-orbit misalignment, rapid rotation, and high j3 may increase 
the spread in debris energy. Finally, we quantify for the first time the gravitational 
wave emission due to the strong compression of a star in a high-/? TDE. Although 
this signal is unlikely to be detectable for disruptions of main sequence stars, the tidal 
disruption of a white dwarf by an intermediate mass black hole can produce a strong 
signal visible to Advanced LIGO at tens of megaparsecs. 



1 INTRODUCTION 



Stars which pass too close to supermassive black holes 
(SMBHs) are disrupted by the enormous gravitational gra- 
dients acting on them. The eventual fallback of ~ 50% of the 
star's mass onto the black hole can produce a highly lumi- 
nous, multiwavelength flare - the primary observable signa- 
ture of a tidal disruption event (TDE). Over a dozen strong 
TDE candidates have been observed, with most detections 
made in X-ray ( Bade et al.|1996 Komossa fc Greiner|1999 



on iVxDE between 10 * and 10 ^ yr ^ ( 


Frank & Rees|l976 


Lightman & Shapiro||1977| |Cohn & Kulsrud |1978| |Magor-| 


rian & Tremaine||1999 


Wang & Merritt||2004|). Alternative 



mechanisms for enhancing the TDE rate include triaxiality 
in a galaxy's nuclear potential ( Merritt fc Poon||2004[ ), en- 
counters with massive perturbers ( [Perets et al.||2007|), the 
effect of an inspiraling secondary SMBH (Ivanov et al. 12005 



but some in optical archival (van Velzen et al 



transient (Cenko et al. 2012 Gezari et al. 2012 



2011[ ) and 
searches. 

Recently, coUimated jets from two relativistic TDE candi- 



Gezari et al.]|2003| ) or UV |Gezari et al.||2006| |2008| |2009| ), 



Chen et al.|2009[|2011||Wegg fc Nate Bode|2011[ ), and gra7 
itational wave recoil after the merger of a binary SMBH 
( Stone fc Loeb|2011[ ). 

Past theoretical work on TDEs has included analytic 
estimates of event energetics and timescales (|Rees| 1988| 



dates have been detected by the Swift satellite ( Levan et al 



2011 



2012 



[Bloom et al.||2011| [Zauderer et al7||20TT| [Ceiiko et al 



Phinney||1989| |Ulmer||1999[ [Strubbe fc Quataert||2009[ ), hy- 
drodynamical simulations of the disruption process in both 



The rate of TDEs is highly uncertain, from both theo- 
retical and observational perspectives. On the observational 
side, uncertainties stem from both the low sample size and 
unclear sources of systematic error; nonetheless, observa- 
tional estimates of the TDE rate per galaxy g enerally flnd 



smoothed-particle- hydrodynamics ( 


Nolthenius fc Katz 1982 


Evans fc Kochanek|| 


19891 ILaguna 


et 


al. 1993 Lodato et 


al. 


2009|) and mesh ( 


Guillochon et al. 


2009 


Guillochon fc 


Ramirez-Ruiz||2012|) codes, and radiative transfer work to 



10" 



10"^ yr~^ (Donley et al. 



2002, Gezari et 



al.|2008 1. This is in rough agreement with the wide range of 



quantify emission and absorption processes in TDEs ( Kasen 
fc Ramirez-Ruiz|2010| [Strubbe fc Quataert|2011| ). The large 
hierarchy of time and length scales involved in a TDE makes 
it difficult to self-consistently simulate one from disrup- 
tion to the onset of accretion, so other work has focused 



theoretical predictions for the TDE rate, which invoke differ- 
ent ways to scatter stars into the SMBH "loss cone" (the re- 
gion of low angular momentum phase space containing orbits 
which pass inside the tidal sphere). The most theoretically 
secure method of feeding stars into the loss cone is standard 



on the formation (Kochanek 1994) and evolution (Mon- 



tesinos Armijo & de Freitas Pacheco 



20111 of TDE accretion 



disks. The properties of these flares depend crucially on the 
spread in orbital speciflc energy of the post-disruption de- 
bris streams, as that spread sets the mass fallback rate onto 
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the SMBH. All fallback rates are generally expected to pro- 
duce multicolor blackbody emission from the accretion disks 
they feed (Lodato & Rossi 20111, while super-Eddington 



dependences on stellar structure (|Diener et al. 19951, stel 



fallback rates may drive powerful outflows, increasing the 



optical luminosity by orders of magnitude ( Loeb & Ulmer 
TI97I IStrubbe & Quataert''2009l [SoTTl). 



Much of the pioneering work on TDEs was done in the 



1980s using the semi-analytic afhne model (Carter & Lu 



minet 



1983 hereafter CL83), which treats the disrupting 



star as a set of concentric ellipsoidal shells evolving under the 
combined influences of self-gravity, pressure, and the SMBH 



tidal field ( Carter & Luminet 1985 Luminet & Carter 1986 1. 



This model has found a wide range of uses, and has been 
generalized to include both thermonuclear reaction networks 



I Luminet fc Pichon||l989 1 and general relativistic (GR) ef- 
fects ( Luminet fc Marck|1985 1, although its validity tends to 
break down at late times as the stellar debris exits the tidal 
sphere. One key finding of the affine model is that during 
the early stages of disruption, prior to the star's arrival at 
pericenter, motion orthogonal to the orbital plane decouples 
from motion within, leading to a strong, one-dimensional 
compression (a vertical "pancaking") of the star. This ef- 
fect scales strongly with the penetration factor /3, defined 
as the ratio of the tidal radius Rt to the pericenter radius 
Rp. The pancaking of the star is reversed by the buildup of 
internal pressure, which leads to a rebound in the vertical 
direction. Shock formation accompanies this rebound (and 
occasionally the infall period prior to maximum compres- 
sion), with X-ray shock breakout a potential though as yet 
undetected observational signature of TDEs ( jKobayashi et| 
al][2004l [Guillochon et al.|2009| ). 

In this paper, after establishing basic dynamical fea- 
tures of TDEs (§2) we present a new analytic model to an- 
alyze the tidal free fall of the star prior to its maximum 
vertical compression (§3). In many ways, this represents a 
simplification of the affine model, and its primary appeal is 
its greater analytic tractability. Using our model we correct 
a longstanding error in the literature on the spread in debris 
energy Ae. We verify the robustness of our estimates by con- 
sidering redistribution of vertical collapse energy to in-plane 
motion (§4), the desynchronization of vertical collapse (§5), 
and leading-order GR corrections (§6), the latter of which 
are found to be small. We examine the gravitational waves 
(GWs) generated by rapid changes in the star's quadrupole 
moment during maximum compression, and find them to 
be detectable by Advanced LIGO for disruptions of white 
dwarfs (§7). We conclude with the observational implica- 
tions of our work, which primarily involve the suppression 
of strongly super-Eddington TDEs (§8), and a general dis- 
cussion (§9). 



2 DYNAMICAL ENERGY SPREAD 

A star is tidally disrupted if the pericenter of its orbit, Rp, 
lies inward of the tidal radius, 



Rt = R,{Mbh/M, 



,1/3 



(1) 



Here M, and R, are the mass and radius of the victim star, 
and Mbh is the black hole mass. In reality, this expression for 
the tidal radius is not exact and contains weak, order unity 



lar spin, and black hole spin ( jKesdenj2011 1. We ignore these 
complications in this paper. Very shortly after entry into the 
tidal sphere (and before pericenter passage for /3 > 1), the 
SMBH's tidal forces do an amount of work exceeding the 
star's gravitational binding energy, and the star's fluid ele- 
ments begin moving on roughly geodesic trajectories. In the 
standard picture, hydrodynamic forces are subsequently ne- 
glected and the specific orbital energy e of the debris streams 
is "frozen in," with a spread given by 



Ae^k 



GMbhR* 



(2) 



ApTp 

Ae = VpAVp 



where G is the gravitational constant and k a constant of 
order unity related to stellar structure and rotation prior 
to disruption. This approximate estimate can be obtained 
by Taylor-expanding the SMBH potential around the star 
at pericenter, or alternatively by multiplying the equivalent 
tidal acceleration at pericenter Ap ~ {GMbu/ R^){R* / Rp) 
by the dynamical time Tp ~ (GMbh/-Rp)"^''^ to get AVp = 
Using Vp = (2GMBH/-Rp)^''^ one can then find 
GMBYi R*/Rl. We note that Eq.j2|is widely 
used in the literature (| Evans fc Kochanek| 1989| [Kochanek 
1994; Ulmcr 1999, Kascn fc Ramircz- Ruiz|2010| [Strubbe fc 
|Quataertt2009| , Lodato fc Rossi,2011J . 

However, this reasoning is incorrect; by the time the 
star reaches pericenter its fiuid elements are moving on al- 
most ballistic trajectories. As the star plunges into the tidal 
sphere, internal forces become subdominant to the SMBH 
tidal field, with the ratio of the tidal to the self-gravitational 
acceleration given by at/og « (Rt/R)^. Here 7? is the orbital 
separation. The work done by internal forces decreases more 
slowly, ~ GMtR/ R^, although this simple expression over- 
estimates the amount of work done by internal forces, which 
at R ~ Rt will self-cancel each other to first order (given 
that the star is initially in hydrostatic equilibrium). 

To accurately evaluate Ae at pericenter passage, one 
would need to account for distortions in the free-falling star's 
physical shape, as well as internal velocities. At any point 



along the star's orbit, a Cartesian coordinate system (Bras- 
sart fc Luminet||2008 \ will define the principal axes (eigen- 



vectors) of the tidal tensor. If we define X parallel to the 
vector connecting the star and SMBH, Y in the orbital plane 
but perpendicular to X, and Z perpendicular to the orbital 
plane, the star will be stretched in the X direction but com- 
pressed in the Y and Z directions. The compression along 
the Y axis reduces the potential gradient across the star, 
invalidating the above formula; further inaccuracy is intro- 
duced by the internal motions (i.e. velocity shear among 
ballistic debris trajectories) of the star within the SMBH's 
tidal sphere. A more accurate estimate of the spread in spe- 
cific energy can be found by taking the potential gradient at 
the moment of tidal disruption, i.e. when the star crosses the 
tidal sphere and becomes unbound, as after this point the 
motion of the debris becomes roughly geodesic. This revision 
to the approximation of energy freeze-in yields 



k 



GMbhR* 



(3) 



We note that an analogous conclusion (on the /3- 
independence of the energetics of tidal disruption) can be 
seen in tidal separations of binary stars by SMBHs (Sari 



El 0000 RAS, MNRAS 000, 000-000 



Consequences of Strong Compression in Tidal Disruption Events 3 



et al. 2010 hereafter SKRIO). We can alternatively use 
(3 = Rt/Rp to rewrite Eqs.Illslas 



k0 



.„ G MbhR* 



(4) 



with n = 2 for the standard, Eq. |2] picture and n = for 
our revised, Eq. [3] analysis. In the sections below, our more 
detailed analysis of the tidal compression experienced by the 
star will examine if intermediate or piecewise values of n are 
more appropriate. The observational implications of changes 
to Ae are discussed in 58. 



3 FREE SOLUTIONS, AND FREE COLLAPSE 

Other factors could influence or invalidate the simple an- 
alytic argument presented in §2, such as redistribution of 
energy during the moments of maximum compression, GR 
corrections, stellar spin, or simply work done on the star's 
fluid elements by subdominant internal forces inside the tidal 
sphere. In this section, we introduce a new analytic model 
for the tidal free fall of a disrupted star that will help us 
approach these issues. 

Because the dominant source of TDEs is expected to be 
stars scattered onto radial orbits from 



pc scales (Magor- 
|rian fc Tremaine||1999l |Wang fc Merritt||2004[ ), we assume 
a parabolic orbit for the center of mass of the star, with 
distance from the SMBH given by: 



R = 



2Rp 



(5) 



1 + cos / 

For such an orbit time t is related to true anomaly / via 

1 / 2R^ \^^^ 
t = I ( gMbh ) ^^""^^/^^ + tan^(//2)) , (6) 



although the differential form 



dt 



81/2 



(1 + cos/)' 



GMbh 



(7) 



/ 



is more generally useful. We set t 
use / > throughout this paper. 

The sequence of events in a TDE, first noted by CL83, 
will be useful shorthand for us, so we introduce it here. Phase 
I (near-equilibrium) of a TDE lasts while R > Rt, and the 
star remains in approximate if slightly perturbed equilib- 
rium. Phase II (free fall) begins when the star crosses the 
tidal sphere and becomes gravitationally unbound; in this 
paper we will treat the transition between Phases I and II 
as instantaneous, an assumption we justify below in §5. The 
assumption of tidal free fall is very useful because of the 
existence of analytic, "free" solutions to the Hill equations 
in the parabolic restricted 3-body problem, but it is not im- 
mediately obvious that internal forces in the star can be ne- 
glected for R < Rt. To first order the approximation seems 
reasonable because the ratio of tidal acceleration to self- 
gravitational acceleration grows quickly, as at/flg ~ (Rt/R)^ 
for the bulk of the star. Furthermore, the star's internal 
pressure and self-gravity partially cancel each other, further 
reducing their combined contribution. For now, we assume 
the validity of the free fall assumption, but after developing 
further machinery we will justify it further in §5. 



During this free fall, the star is compressed perpendicu- 
lar to the orbital plane (along z) and in one direction within 
the orbital plane, while being stretched along the other in- 
plane direction. Although for the limiting case of radial infall 
the problem is self-similar in all three dimensions (SKRIO), 
the rotation of the line connecting the star's center of mass 
to the SMBH breaks the in-plane similarity. By the time the 
star has reached pericenter, the x direction (which is parallel 
to the line between the SMBH and the orbital pericenter), 
is compressed, and y is stretched, but the distortions are 
both much smaller than the compression orthogonal to the 
orbital plane. Shortly after passing pericenter, synchronous 
tidal free fall in the z direction leads to very strong com- 
pression of the star, which is eventually reversed by hydro- 
dynamic forces. Phase III (bounce) begins when hydrody- 
namical forces become strong enough to begin slowing the 
collapse of the star along its vertical axis. Once the star's 
collapse has reversed, hydrodynamical forces quickly become 
negligible again, and Phase IV (the rebound) begins, with 
stellar gas once again moving on ballistic trajectories. 

We take as initial conditions for Phase II a spherically 
symmetric star at the tidal sphere, with fluid elements pos- 
sessing initial positions r (the coordinate origin tracks the 
star's center of mass) and initial velocities in the center of 
mass frame u{r). Making the approximation that upon en- 
tering the tidal sphere, internal forces become negligible un- 
less and until compression triggers shock formation or isen- 
tropic pressure buildup, we take the pre-shock trajectories 
of these fluid elements to be completely ballistic. This means 
that their trajectories are given by the parabolic Hill equa- 
tions, fu = {xn,yH, zu}- The free solutions to these equa- 
tions, neglecting self-gravity, can be written in closed form 
(SKRIO) using coordinates where distance has been normal- 
ized by _Rt and time by \jRi/ (GM,)\ we denote such coor- 
dinates in this paper by writing tildes over them. All other 
coordinates are in physical units, unless otherwise noted. Al- 
though there are 6 independent solutions to these equations, 
motion out of the orbital plane is decoupled from motion 
within it, so only two are relevant for perturbed motion in 
the z direction: 

5h — EzE + FzF 
1 2 sin / 



Ze 



Zf = 



/3 1 cos / 
1 2 cos / 



(8) 



/3 1 + cos / 

Here E and F are undetermined coefficients that are set by 
the initial conditions described above. In particular, if we 
require that a fiuid element of initial position z = zo has 
initial velocity u' = i = Oat/ = /t, where the true anomaly 
upon entry into the tidal sphere is given by 



ft 



■ arccos(2//3 — 1), 



then 



E - 
F ■- 



'Zo 



zo. 



(9) 

(10) 
(11) 



If we introduce a tidal potential felt in the rest frame of 
the star, the tidal acceleration is given by SKRIO as 



dz 



_^3(1+C0S/)!-, 



(12) 
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We note that the self-similarity of Eq. [12] implies that the 
free solutions all collapse to z — simultaneously at a true 
anomaly fc, although physically this collapse will be re- 
versed shortly before by the buildup of pressure gradients 
strong enough to counteract the tidal forces compressing 
the star. However, it is useful to solve for /c using Eq. [8] 

From this formula we see that in the limit of /3 — ^ oo, col- 
lapse along the z-axis occurs at fc — 0, i.e. at pericenter, 
while in the marginal disruption limit of /3 — > 1, collapse 
occurs at fc — 7i"/2, i.e. at a fixed point past pericenter. We 
have already mentioned that the free solutions become less 
valid for small /? due to the increasing importance of inter- 
nal forces, but we can see from Eq. [13] a second, stronger, 
inconsistency at low /?, which is that the free solutions dic- 
tate vertical collapse after the disrupted star leaves the tidal 
sphere, i.e. fc > \ft\- This occurs for /3 < 1.3. 

Although the onset of Phase III is dictated by compres- 
sion in the z direction, the outcome of the bounce will be 
affected by motion within the orbital plane during Phase II, 
when f < fc- We therefore describe here the free solutions 
within the orbital plane (SKRIO): 

in = Ax A + Bxb + Cxc + Dxd 
yu = AyA + ByB + Cyc + Dyu 
1 sin / 
" /3 1 + cos / 
1 cos / 



XA 



VA = 



Xb 



/3 1 + cos / 

1 ■ t 

= sm r 

/3 ^ 
(1 + cos/) 



(14) 



cos/) 



yc = 



/3 



cos / tan(//2) 



X-D 



(8 + 12 cos/) tan (//2) 



2/D = 



1 35 sin / - 2 sin(2/) -h 3 sin(3/) 



P (l+C0s/)2 

If we consider a point on the star with an initial po- 
sition, relative to the star's center of mass, of (xo, 2/0,20) 
and zero initial velocity (here, as before, "initial" refers to 
f = ft, i.e. crossing into the tidal sphere), then we have 4 
initial conditions for 4 unknowns: {A, B, C, D} . Using Eq. 
[9] we find 



A : 

B 

C 
D 



= ^ (-8iov/^^ + 2yo(/3' + 2/?-4)) 
= ^(2iov/^(/3'- 4/3^ + 8) 

-hyo(9/3^- 12/3^-8/3 + 16)) 
= ^ (io(2/3' +/3-2)- 2yo^P^{p' - 1)) 



(15) 
(16) 

(17) 
(18) 



l.Orr 

0.5- 
^ 0.0- 
-0.5- 
-1.0 + 
0.5- 
« 0.0- 
-0.5- 
-l.OU 




(a) 



(c) 



-T ■■ iV 



(b) 



I i . i i i n i i i I i i i i I i i i i I i i i i I i i i i I 



(d) 



/ 



Figure 1. Normalized height z = z/i?* versus true anomaly / 
for the vertical collapse of one-dimensional stars at varying /3. 
The solid blue curves in panel (a) are /3 = 2; the dashed orange 
curves in panel (b) are /3 = 4; the dotted purple curves in panel 
(c) are /3 = 10; the dot-dashed brown curves in panel (d) are 
/3 = 40. Each scenario is initialized at / = /t(/3). Note that / = 
corresponds to pericenter. 



All six of the free solutions we have listed can be thought 
of as slight perturbations to different orbital elements of the 
parabolic center of mass trajectory, boosted into the center 
of mass frame. 

We have now exactly specified the motion of the ideal- 
ized star's fluid elements in the orbital plane during Phase 
II. We plot the vertical free solutions for a variety of /3 in 
Fig. [1] and snapshots from motion within the orbital plane 
in Fig. [2] Here we list several important features of the free 
solutions, when they are initialized with static spheres of 
matter at / — ft: 

• For / > /t, an initially spherical shell of matter will 
deform into a sequence of roughly ellipsoidal shapes. It is 
simple to demonstrate that they do not generally take the 
form of true ellipsoids, however. 

• Initially concentric spherical shells of matter remain 
concentric, with in-plane principal axes that remain aligned 
with those of other concentric shells. 

• Slices of the star through the orbital plane (z=0) main- 
tain reflection symmetry across their rotating in-plane prin- 
cipal axes. 

• The derivation of the free solutions assumes that 
R,/R <C 1 (SKRIO). If we neglect stretching of the star, 
this is equivalent to requiring /3 <^ (A/*/A/bh)^'''', a condi- 
tion that is in general easily satisfied: a W^Mq SMBH, if 
non-spinning, cannot disrupt solar- type stars with /3 > 11 
(they will plunge directly into the horizon). Even a maxi- 
mally spinning SMBH of this mass cannot disrupt solar-type 
stars with j3 > 47. The effects of tidal stretching will make 
it somewhat harder to satisfy this assumption, but only for 
the minority of the star's mass that is strongly stretched. 

The free solutions allow us to directly solve for the stel- 
lar axis ratio as a function of / or t, and it is trivial to do so 
numerically, but there is an exact analytic solution as well. If 
we denote the lengths of the long and short principal axes of 
our tidally distorted star (within the orbital plane) as riong 
and Tshort, respectively, we can solve for them by rewrit- 
ing xo — cos 60, yo = sin 60, and finding the appropriate 60 ■ 
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Figure 2. The x and y coordinates (in units of Rt, which for the IO^Mq SMBH in this example is 100 Rq) of the free solutions for 
varying /3, translated so that the origin lies on the SMBH. As before, we mark the /3 = 2 trajectory as solid blue, /3 = 4 as dashed orange, 
/3 = 10 as dotted purple, and /3 = 40 as dot-dashed brown. The free solutions for an initially circular midplane slice of a star are magenta 
at / = ft, pink at / = 0.7/t, and red at / = fc- The tidal radius is marked as a gray dashed circle. The right plot is a zoomed-in version 
of the left. The free solutions are breaking down for the /9 = 40 curve near pericenter, as the long axis of the star exceeds the orbital 
radius in size. 



More specifically, we set -^Ruif) = (with = a^H + J/n)) 
and solve for ^cx, the values of which extremize -Rh- More 
physically, we are searching for the initial angles 6'ox around 
the star which at a later orbital phase / will correspond to 
its principal axes in the orbital plane. Once we have the ini- 
tial angular positions of the principal axes, Sex, we can plug 



in to Eqs. (14 1 and solve for the size of the principal axes 
at a later true anomaly f > ft- We also find the misalign- 
ment angle f between the long in-plane principal axis and 
the orbital velocity vector. The in-plane stellar geometry is 
presented in Fig. |3] 

The algebra involved in this solution is unenlightening, 
so we leave the general solution Sex(/) for numerical work 
and only derive analytic expressions for Sex(/c), which is 
the situation of greatest interest. The details are contained 
in Appendix A, but we plot the results below in Fig. [4] 
Specifically, these are the sizes of the principal axes at / = 
/c. For comparison we plot curves of the high /3 limiting 
behavior, for which fiong ~ ^P^^^ + ^/3~^^^ and f short ~ 



2/3- 



-1/2 , 



-3/2 



The primary interesting feature of the axis ratio calcu- 
lations is that for disruptions of stars by supermassive black 
holes, the physically relevant range of riong and rshort is quite 
narrow, being confined between 3 to 5 for the former, and 
0.3 to 0.5 for the latter. For the tidal disruption of a star by 
an intermediate mass black hole (IMBH), or perhaps more 
exotic disruption scenarios, a larger range of P (and there- 
fore riong, ^short) cau be attained, but for star-SMBH TDEs 
only a surprisingly narrow range of principal axis lengths are 



accessible. This implies that the naive Taylor expansion of 
the SMBH potential as the star passes through R — Rp, i.e. 
Eq. [2j will fail primarily because of internal velocities within 
the free-falling stellar debris, and only secondarily because 
of distortions in the star's shape. 

As we shall see in §4, when estimating energy redistri- 
bution during maximum vertical compression, the misalign- 
ment angle i/ plays a larger role than the slowly-varying 
axis ratio. This will prove relevant when calculating cor- 
rections to Ae, and is shown in Fig. [s] The angle Vc (as 
elsewhere, the subscript c denotes evaluation at / = fc) is 
found to be a rapidly decreasing function of /3; to leading 
order, tani^c oc I3~^^^. 

We can now use the free solutions {A, B, C, D, E, F} to 
quantify precisely the spread in debris energy at the tidal ra- 
dius. Because these solutions can be thought of physically as 
perturbations to the orbital elements of a parabolic trajec- 
tory, all possess exactly zero energy except for the in-plane 
"D" solution, which has specific energy given by 



20GM.D /Mbh 



R, 



M, 



1/3 



/3, 



(19) 



where D is the coefficient of the fourth in-plane free solu- 
tion, corresponding to slight variations in the eccentricity 
of a near-parabolic orbit (SKRIO). If we initialize our free 
solutions with an unperturbed sphere, i.e. Eqs. |15| then we 
find a specific energy for each fiuid element of 

e„ = g^feH-^' (^o(i _ 2/13) + 2yoVl3-' - P'^), (20) 
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-4 -2 2 4 

X 

Figure 3. An initially (/ = ft) circular ring of stellar fluid ele- 
ments has been tidally distorted by the time it reaches f = fc- 
The principal axes of the distorted, free-falling body are the solid 
green vectors, the center of mass velocity is the dashed yellow 
vector, and the direction to the SMBH is the dotted red vector. 
The angle tpc (Tc) is measured between the negative x direction 
and the long principal axis (stellar velocity vector). We define the 
misalignment angle Uc = Tc — V'c 




Figure 4. The principal axis lengths, Hong and rghort, of the 
distorted star (at / = /c) vs /3. Here rj^ng is plotted as thick green 
curves; rghort thin yellow curves. The exact solutions are solid 
lines, while the dashed curves are the high-/3 Taylor expansions 
given by Eqs. ( |A6| , \A7\ . 



where xo and yo are the initial positions of a debris stream 
relative to the star's center of mass at ii = _Rt normalized 
by the stellar radius. Notably, the specific energy is to lead- 
ing order independent of /3, with the weak /3-dependence 
becoming negligible at high penetration factors. Defining 
Xo = fo cos (po and yo — fo sin 00 , we can analytically ex- 
tremize Eq. [20] with respect to </io, to find that the spread in 
energy of these unperturbed free solutions is actually fully 




Figure 5. Curves illustrating misalignment of the tidally raised 
bulge and the orbital velocity vector at the point of z-coUapse 
(/ = fc)- The red dotted curve is the angle between the negative 
X-axis and the tidal bulge (V'c), the green dashed curve is the 
angle between the negative x-axis and the orbital velocity vector, 
and the solid yellow curve is the difference between them, i.e. 
the misalignment angle (uc)- These angles are plotted against the 
penetration factor /3. The high /3 limit for Uc is the dot-dashed 
yellow curve. 



independent of /3: 



(21) 



In our idealized model of a spherical, stationary star under- 
going tidal free fall, the assumption of energy freeze-in at 
the moment of disruption implies n = 0, and, surprisingly, 
not even the weak /3 dependence one might expect from Eq. 

Eol 

In the following three sections, we examine the robust- 
ness of this model, and consider possible corrections to our 
expressions for Ae. With limited exceptions, we find that 
the arguments made in this section remain generally valid. 



4 TOTAL VERTICAL COLLAPSE AND 
BOUNCE 

For f ~ fc, motion in the vertical direction has decoupled 
from in-plane motion and the star undergoes a homologous 
vertical collapse. In this regime, the vertical velocity of the 
free solutions near the point of maximum collapse is very 
close to a constant value, with Wc oc /3, a result known since 
CL83, although the exact value, for arbitrary fc, is 



1/2 



(l_/3-l)V2_^l 



(22) 



With this formula we can begin thinking about Phase III 
of a tidal disruption event, and in particular whether it can 
alter Eq. [2l] 

As we have seen in the previous section, once f ~ fc, 
the majority of the star simultaneously "pancakes" into a 
sheet of matter strongly compressed in the vertical direc- 
tion. If non-gravitational forces were truly negligible, an ide- 
alized, one-dimensional {z extent only) star would momen- 
tarily possess zero height at f = fc, but in reality sufficient 
compression will create a pressure gradient strong enough to 
oppose free fall in the z-direction. The resulting bounce will 
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reverse the vertical free fall and lead to vertical expansion at 
speeds comparable to Wc- The vertical rebound will have a 
limited impact on Ae because it is effectively decelerated by 
the tidal potential (for example, if we generously approxi- 
mate the rebound as elastic due by reflecting w at / = /c, the 
asymptotic free solution velocity w — > as / — > tt), but the 
smaller rebound velocities Av^, Avy in the orbital plane can 
in principle have more significant effects, as Ae ~ VpAv. In 
this section we assume the bounce is adiabatic; in particular, 
we neglect both dissipation in shocks and the thermonuclear 
energy release from the compression of the stellar core. For a 
more thorough discussion of these possibilities see |Luminet| 
& Pichon||T^9|;|Brassart & Luminetl (l2008l. 



During Phase III of a TDE, the requirement that central 
pressure rises to halt the kinetic energy of collapse implies 
that the star's peak internal specific energy will be 



t/c~/3'l7.(\/l-/9-i + l)' 



(23) 



where C/* = GMt/R*. Assuming a polytropic equation of 
state P = Kp'^ , and furthermore that strong compression 
in the z direction means that the density enhancement will 
be due to collapse in z alone (since the cross-sectional area 
within the orbital plane, « vrriongrshort remains roughly con- 
stant), gives a minimum stellar height and maximum stellar 
density of 



Zmin 

1?7 



Pc 



2/(7-1) 



(24) 



where p« is the mean pre-disruption stellar density. The du- 
ration of maximum compression is a steep power of the im- 
pact parameter, specifically 



^13 



-(7 + l)/(7-l) 



r., 



(25) 



with r, = 1/^Gp*. Although we have only derived these 
formulas at the order of magnitude level, they have been cal- 
ibrated over a wide range of /? by both the affine model ( |Lu-| 
minet fc Carter|1986 l and one-dimensional hydrodynamical 
simulations (Brassart &: Luminet 20081 (hereafter BL08). 



Specifically, for 7 = 5/3 polytropes, the affine model found 
[/c « 1.2^7.;g^ Pc = 1.3p./3^ and = 8.5r.^~*, calibra- 
tions which were essentially duplicated in BL08. Likewise, 
Zmin ~ 4.5^-^i?, if we assume that the rise in density comes 
entirely from homologous, vertical stellar collapse (i.e. that 
the in-plane area of the compressed star is ~ riongrshort). 

If we assume that the pressure-driven bounce acts 
isotropically (i.e. that shear stresses from viscosity or shocks 
remain unimportant), then the relevant changes in velocity 
can be estimated as pcAvi/rc ~ APc/n, with n the physi- 
cal dimension of the star parallel to which pressure gradients 
impart Avi. Specifically, 



Av^ 

AWshort • Vc 
A«lonK • Vc 



'Uc 



/Uc sin(i/c) 

'^short 



^Uc COS(l/c), 



(26) 



ong 



where we have denoted the direction of the orbital velocity 
at / = /c as the dimensionless unit vector Vc- 

This leads to energy perturbations within the orbital 
plane at bounce of Aem ~ VcAv. Using our exact formulae 
for axis lengths and alignment, we plot the results in Fig. 




Figure 6. Fractional specific energy perturbations during the 
bounce phase. We plot perturbations along the short axis of the 
star as thick yellow curves, and along the short axis as thin green 
curves. Exact solutions are solid lines, and the leading order be- 
havior from Eqs. \27\ and | |28[ l are dashed lines. The bounce rep- 
resents at most a factor 2 correction to specific energy of the 
stellar debris for 4 > /3 > 3; above these values, the bounce is 
negligible. Below /3 3, a larger correction is possible, but the 
free solutions become somewhat unreliable. The high degree of 
alignment between the stellar bulge and the orbital velocity vec- 
tor causes perturbations along the long axis to actually dominate 
those along the short axis above /3 15. 



[6] along with the limiting behavior at high /?, which is well- 
approximated by the Taylor expansions in Appendix A as: 



Aeiii,] 



Aeiii , 



1 + 



' 31^ 
II 



-5/2 



Aeu 



Aeu. 



(27) 
(28) 



Here we have assumed 7 = 5/3, and that Zmin ~ 57?«/3~'^ 
based on the BL08 calibration. We have also approximated 
Vc ~Vp, which is accurate for high /3 though a mild overesti- 
mate at low p. Even with this overestimate, we can see from 
Fig.[6]that only for /3 < 3 (where our model's assumption of 
tidal free fall begins to break down) can the pressure-driven 
bounce along the short principal axis of the star provide an 
order unity enhancement to the total spread in (in-plane) de- 
bris energy. The contribution of the bounce along the longer 
principal axis remains negligible at all (3. 

Because the z-component of 14 = 0, the spread in ki- 
netic energy of vertical motion at the time of bounce is given 
by 



Aeiii , 



Avt 



(29) 



Interestingly, for P > 10 the total spread in kinetic energy at 
the time of bounce is dominated by Aem.z, not Aeu. How- 
ever, the instantaneous vertical kinetic energy at f — fc 
will disappear as / — >■ vr. Even neglecting dissipation of the 
kinetic energy of vertical free fall into shocks, and assum- 
ing a perfectly elastic bounce, the tidal potential (Eq. 12 1 



will efficiently decelerate the vertical motion of the debris 
during the phase IV rebound and later expansion. This can 
be seen by continuing the free solutions past f ~ fc, which 
corresponds to a reflection of vertical velocity. 

Here we have ignored energy release from thermonu- 
clear burning at the time of maximum compression, which in 
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principle could increase Ae. However, past estimates made 



in the framework of the afhne model (Luminet & Pichon 



1989 Table 8) found that for the range of /3 values consid- 
ered (5 ^ /3 20), the total thermonuclear energy release 
was less than Uc, making it unlikely to change the analysis 
of this section. 

Of course, our analysis of energy redistribution in Phase 
III depends critically on how synchronously the vertical col- 
lapse of the star proceeds. If individual "columns" of the star 
do not collapse in the homologous manner implied by Eq. |12[ 
it is unlikely Zmin will reach the extreme values predicted by 
the simple arguments in this section. Alternatively, if sepa- 
rate columns collapse in a desynchronized way, it is possible 
that pressure waves from collapsed regions of the star will 
propagate upstream to uncoUapsed regions and cause them 
to rebound prematurely. In either scenario, the effective Zmin 
will be enhanced, enabling greater coupling of the bounce en- 
ergy to motions within the orbital plane, and increasing the 
values of Aein, short and Aemjong. Therefore, Eg. |21| should 
be regarded as a lower bound on Ae - a higher value of n 
would be favored if the desynchronization of vertical collapse 
transfers kinetic energy to in-plane motions more efficiently 
than in our estimates here. A similar effect should arise in 
hydrodynamical simulations of TDEs that lack sufficient ver- 
tical resolution to capture the maximum compression of the 
star (we discuss this further in §9). In the following section, 
we consider physical sources of desynchronization. 



5 DESYNCHRONIZATION 

The synchronous vertical collapse of a one-dimensional star 
into a thin, pancake-like sheet only occurs if the initial distri- 
bution of vertical velocities is self-similar, i.e. wo{zo) oc zo. In 
previous sections we have assumed the trivial self-similarity 
of Wo = 0. Deviations from self-similarity will be seeded at 
early times by the nonlinear hydrodynamics of actual disrup- 
tion at the tidal radius Rt, and also later, as the self-gravity 
and pressure of the stellar debris perturbs the free solutions 
for / < fc. In this section we quantify in an approximate 
way the efi^ect of desynchronization on our idealized earlier 
conclusions, finding that both the stellar properties during 
Phase III (important for any shock breakout signal) as well 
as Ae could be significantly altered. However, we then argue 
that past hydrodynamical simulations indicate that desyn- 
chronization is likely to be suppressed in physical TDEs, 
justifying our use of the parabolic free solutions. Finally, we 
consider the desynchronization of stellar collapse in three 
dimensions. 



5.1 Desynchronized Free Solutions 

At / = /t, during the transition from Phase I to Phase 
II, velocity perturbations of order ~ ^/G/M^^JR^ could be 
imprinted on the free-falling stellar debris. Normalizing our 
initial conditions {zo, uiq} in units of _R* and GKh / , we 
derive coefficients for the "perturbed" (i.e. wq 7^ 0) vertical 




Figure 7. Here we plot several sets of free solutions (2 versus 
t) for /3 = 2. In panel a, the initial vertical velocity wq = 0. In 
panel b, all fluid elements in the star receive initial velocity per- 
turbations |uio| = \/ GKI, I i?» ; in panel c, |uio| = \J GM, / R:^ /Z. 
In panel d, the star receives homologous velocity perturbations 
uiQ = —zq\J GMt I Rt, causing synchronous collapse before peri- 
center passage. 



free solutions to be 



Ep = -^o^/P^-woJ^{P-2) (30) 



-Fp = 5o + i5o y ^ \/~i3~- 



(31) 



Therefore the true anomaly of a perturbed vertical collapse 
to 2 = is 



tan(/c) 



^o + wo^/27?^/^/^ 



(32) 



We note that both fc and go oc I3~^^^ in the large /3 
limit. Unless Wq oc zq, the collapse will be non-homologous, 
with fc depending on zq. Modest deviations from homolo- 
gous initial conditions will desynchronize the collapse, which 
we illustrate by plotting the desynchronized free solutions 
for 13 — 2 and /3 = 10. We can see that in both cases, the 
time at which the free solutions cross the orbital plane be- 
comes strongly desynchronized, which complicates our pre- 
viously simple treatment of the transition from "tidal free 
fall" to "pressure-driven bounce" and also raises the possi- 
bility that the vertical kinetic energy of free fall could be 
effectively isotropized and transferred to motions within the 
orbital plane, restoring a /3 dependence to Ae. 

We can estimate the amount of desynchronization by 
using trigonometric identities and Eqs. [13] and [32] to find 
A/ = /c-/c. Specifically, 



tan(A/) = 



Ao 



2/3 + \o^/W- 



(33) 



where Xo — wo/ zo . 

Interestingly, both Figs. [7] and [S] show that most of the 
star's desynchronized free solutions have two crossings of 
the orbital plane, raising the possibility of a double bounce 
in desynchronized collapse scenarios (something previously 
seen only due to GR effects, e.g. Luminet & Marck ( 1985 1 



- see §VI). But is it realistic to expect desynchronized col- 
lapse? 



© 0000 RAS, MNRAS 000, 000-000 



Consequences of Strong Compression in Tidal Disruption Events 9 




tidal free fall (BL08, Eq. 43). We note that this is a geo- 
metric proxy for Zmin- 

• The bounce of the collapsing star occurs after pericen- 
ter passage (BL08, Table 5). As shown above, this places 
a strong constraint on the initial velocity perturbations. 
In particular, let us consider a perfectly homologous col- 



Figure 8. The same as Fig. [7] but for /3 = 10. Desynchronization 
is less severe at higher /3. 



lapse for the sake of argument, with wq = — AoZo. Eq. 32 
will only be positive (i.e. bounce after pericenter passage) if 
Ao < (\/2/3\//3 — 1)"'^, a rather small perturbation (A < 0.04 
for /3 = 7, as is relevant here). 

These numerical results indicate that realistic one- 
dimensional stars behave during Phase II much like the 
unperturbed free solutions we presented in §3, supporting 
our earlier assumption that debris energy "freezes in" from 
f — ft down to the bounce, at / = /c. 



5.2 Validity of Free Solutions in One Dimension 

From the above discussion, it is clear that only modest de- 
viations from self-similarity in the initial velocity perturba- 
tions Wo will produce a strongly non- homologous vertical 
collapse at most realistic /3. We can quantify the magnitude 
of the initial velocity perturbations Ao required to signif- 
icantly desynchronize one-dimensional collapse by making 
the approximation (valid for small A/) that the desynchro- 
nization timescale AtiD ~ Af ^R'^/{GMbii)- If we then re- 
quire Atiu < Tc = XcT*, then for a 7 = 5/3 polytrope (with 
Xc ~ 8.5) and inwardly directed velocity perturbations we 
find the condition that 



|Ao| < 



2/3tan(2xc/3-^'"') 



1 y?^tan(2xc/3-5/2) • 



(34) 



The factor of ~ 2 inside the argument of the tangent 
comes from the difference between r* = ^fxjGpl and 
^/WJ{GM7). This condition grows more restrictive as /3 
increases, with the right hand side of Eq. |34] roughly pro- 
portional to /3~^. 

For one-dimensional stellar collapse, high-resolution hy- 
drodynamical simulations indicate that a highly homologous 
collapse is physically realized (BL08). As noted before, this 
is likely due to a combination of two factors: the partial can- 
cellation of stellar pressure with self-gravity, and also that 
at/flg ~ {Rt/Rf. This explanation is supported by past 
investigations of stellar tidal disruption in the afiine ellip- 
soids approximation: for example. Fig. 4 in CL83 shows the 
first order cancellation of pressure and self-gravity for early 
parts of Phase II. Three results of BL08 further support the 
validity of the unperturbed free solutions in Phase II of a 
TDE: 

• The actual collapse of the star is visually homologous 
during Phase II, as seen by the near-linearity of a vertical 
velocity versus height plot at different times (BL08, Fig. 
2). Although the figure deviates slightly from homologous 
collapse at large radii, possibly due to the fact that the low- 
density outermost regions of the star are disrupted slightly 
before the higher density inner regions (like the peeling of 
onion shells), these outer deviations do not appear to affect 
the key dynamics of the bounce. 



The maximum central compression pn 







2/(7-1), 



5.3 Validity of Free Solutions in Three 
Dimensions 

The full problem of tidal disruption is three dimensional, 
and some three dimensional simulations ( Laguna et al.|1993| 



[Guillochon et al.|2009| ) have indicated that one-dimensional 



descriptions of the bounce phase ( Luminet fc Carter||1986 1 
may strongly overestimate the degree of compression. How- 
ever, lack of vertical resolution in the three dimensional 
simulations makes it difficult to interpret the discrepency, 
and some high resolution simulations ( Rosswog et al.|["2009 l 
do find degrees of compression closer to our analytic ex- 
pectations in §4. Although the impact of higher dimen- 
sional effects on Phase III of a TDE will only be resolved 
through higher resolution hydrodynamical simulations, we 
present here a simple analytic argument suggesting that a 
star made of many columns, each undergoing homologous 
collapse, should attain Zmin comparable to one-dimensional 
predictions. 

Three dimensional desynchronization is an important 
effect that cannot be ignored at high /3: the bounce timescale 
Tc ~ 8.5r,/3~^ for 7 — 5/3, while the time it takes the 
bulk of the star to pass across the tidal radius is Atzu ~ 
1.4r,(M,/AfBH)^''^. This implies that for /3 larger than a 
critical value, 



/3d 



1/12 



(35) 



Tc <C AfsD, and the leading edge of the star will collapse 
and rebound well before the trailing edge. If we assume that 
the star is truly in tidal free fall during Phase II, and seed 
velocity perturbations remain as small in three dimensions 
as has been indicated in one dimensional simulations, then 
each column of the star will reach its maximum compression 
at different t but the same /c. 

This /c remains fixed in space, much like a nozzle, as 
the star passes through it. An example of this "tidal noz- 
zle" has been seen in hydrodynamical simulations of tidal 
disruptions of white dwarfs; for example. Fig. 6 of |Rosswog| 
et al. (20091. Even if we assume the maximum compres- 



in accordance with the assumption of fully synchronized 



sion predicted by one-dimensional models of stellar collapse 
(zmin ~ 4.5i?,/3~^, for 7 = 5/3), the sound speed in the stel- 
lar midplane, Cs,c, will remain a small fraction of the stel- 
lar orbital velocity. Specifically Cs,c/Vp ^ /3^/^(M./AfBH)^''^, 
indicating that unphysically large P values are required for 
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pressure waves from the region of maximum compression 
to communicate upstream to the Phase II material. Unless 
three dimensional effects influence earlier stages of a TDE 
(by seeding large perturbations during the transition from 
Phase I to Phase II), it seems unlikely that the star will be 
prevented from reaching the strong compressions suggested 
by models of one-dimensional collapse. Among other things, 
this highlights the importance of thermonuclear network cal- 
culations for high-/3 TDEs ( [Luminet fc Pichon|198"9 1. 

With these caveats in mind, we generalize the work of 
§111 to perturbed in-plane free solutions, i.e. where every 
fluid element at / = /t has initial positions {xo,yo} but also 
initial velocities {uo,vo}. The in-plane coefficients for the 
corresponding "perturbed" free solutions are 



/3 



^ ( - 8ioy^^+ 2yo(/?' + 2/3-4)) 



+ 



2\/2 /-, 



^3/2 (mo (2 - 3/3) + io^P^W - 2)) 
^ (2xo^P^{l3^ - 4/3^ + 8) + yo(9/3' 
-12/3^-8/3 + 16))4-^^(So(/3^-8 
+ 28/3 - 16) -I- 2{;o^//3~ 1(3/3'' - 6/3 -f 8)) 
^(io(2/32 + ,5 - 2) - 2yov/^(/3' - 1) 



(36) 



(37) 



(38) 



^{iiod - 2/3)v//3- l + So(/3 - 1)')) 
^ (io (/3 - 2) -f 2yo (39) 

If we now calculate the perturbed specific energy of the 
free solutions, we find 

(ao(l - 2/13) + 2yo V/3-^ - /3-^ (40) 

- \/27^(So\//3-l + i3o) 

where the initial velocities have been normalized by 
\/(Th'L,JlL,. Again, there is no leading order /3 dependence 
in the specific energy, although a calculation of Aep does 
not find it completely /3-independent as in Eq. |21[ Nonethe- 
less, the assumption of tidal free fall during Phase II im- 
plies clearly that the frozen-in Ae should be, to leading or- 
der, independent of /3. As a simple test case, we now apply 
these perturbed free solutions to a uniformly spinning star, 
with normalized angular velocity uj = lo l^jGM^jRl such 
that (D = 1 is approximately the breakup frequency (and 
spin parallel to orbital angular momentum). In Eq. |40[ we 
relabel io ~ focos4>o, yo = fosin^o, uo = — tlifo sin^o, 
-Do = a)focos(/>o, and then extremize Ae with respect to 0o. 
Results are plotted in Fig. |9] in general, pre-disruption stel- 
lar spin will enhance the energy spread Ae by a small factor, 
< 2. Low P and high u will maximize the energy spread. 

We note here that large stellar spins misaligned with 
the orbital angular momentum vector could have a much 
greater impact on Ae by inducing vertical desynchroniza- 
tion. A thorough investigation of misaligned spin effects is 
beyond the scope of this work, but as an idealized limit- 
ing case we apply Eq. [34] to approximate, as a function of 




Figure 9. The enhancement Aep/Acu to the energy spread for 
initially unperturbed free solutions, when pre-disruption stellar 
spin (along an axis parallel to orbital angular momentum) is 
considered. The energy spread is plotted against pre-disruption 
spin ui , where Co is stellar spin normalized by the breakup spin 
in previous plots, the solid blue, dashed orange, 
dotted purple, and dot-dashed cyan curves represent /3 = 2, /3 = 
4, /3 = 10, and /3 = 40, respectively. 



/3, the maximum stellar spin allowed before the Phase III 
bounce would be vertically desynchronized. Specifically, we 
set itJo = 1^0 '"o- We plot these results in Fig. |10| and find 
that combinations of high /3 and relatively rapid values of 
stellar spin are needed to strongly desynchronize the vertical 
collapse. 

Fig. [9] can be taken as representative of the effects of 
both primordial stellar spin, and the angular momentum 
imparted during tidal spin-up of the star prior to its full 
disruption, during the transition between phases I and II of 
a TDE. Tidal spin-up is unlikely to produce misaligned spin, 
however, so the more dramatic type of desynchronization 
suggested by Fig. [lO] can only come from the star's original, 
pre-disruption spin. 

For one dimensional stellar collapse, the frozen-in Ae 
will dominate the post-bounce Ae for all /3. For three dimen- 
sional collapse, the numerical literature is less clear, but we 
have argued here that three-dimensional effects are unlikely 
to strongly redistribute energy to in-plane motion, with the 
possible exception of when sufficiently rapid stellar spin is 
misaligned with the orbital plane, or perhaps when a sim- 
ilar misalignment between SMBH spin and orbital angular 
momentum exists. 



6 GENERAL RELATIVISTIC CORRECTIONS 

The results of all prior sections have assumed purely New- 
tonian gravity; however, tidal disruption occurs at an or- 
bital distance Rt < 50i?g, where the gravitational radius 
7?g = GMsYi/ (? ■ At these small separations, ballistic mo- 
tion follows the geodesies of the Schwarzschild or Kerr metric 
rather than free-fall trajectories in Newtonian gravity: gen- 
eral relativity (GR) is important. A fully relativistic anal- 
ysis of the problem of tidal disruption is beyond the scope 
of this paper, although it has been treated in the past in 
the case of the afBne model ( Luminet &: Marckjl985J, and 
in one-dimensional hydrodynamical simulations ([Brassart &: 
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Figure 10. The maximum value of normalized stellar spin uj that 
will not produce significant ID desynchronization leading into the 
phase III bounce. We plot the exact value calculated from Eg. |34| 
as a solid green line, and the asymptotic behavior u) < 2\/2Xc/3~'^ 
as a dashed yellow line. Regions above the curves will experience 
desynchronization of vertical collapse. 



Luminet|2010 1. If we treat the internal dynamics of the star 
in a Newtonian way (i.e. assume tangentially flat space-time 
in the small region occupied by the star), then Eq. 



12 



s de- 



scription of vertical collapse will be modified, to become 



dz 



1 + 3 



(41) 



where L is the orbital angular momentum of the star, R is 
the orbital radius of the star (both in geometrized units) and 
we have limited ourselves to non-spinning black holes. The 



qualitative results of both Luminet & Marck (19851; Bras- 



sart & Luminet (20101 were that the increased strength of 



the GR tidal field (relative to Newtonian gravity) can ac- 
tually result in multiple vertical collapses, each followed by 
separate bounces which are reversed by the relativistically 
enhanced tidal field. For all but the most deeply-plunging 
TDEs {P > 30), the maximum compression is obtained 
on the first vertical collapse and is similar to the Newto- 
nian Zmin ( Luminet fc Marck|1985 Fig. 10). Therefore, even 
though the formation of multiple outgoing Shockwaves could 
be an important outcome of relativistic compression, the 
first-order spread in debris energy is unlikely to be affected 
by multiple compressions for /3 < 30. 

A separate relativistic effect concerns modifications to 
the pre-bounce spread in debris energy, Ae. Eq. [3] was de- 
rived by Taylor expanding the Newtonian gravitational po- 
tential about the star's position when it crossed into the tidal 
sphere, then subtracting the zeroth-order component. We 
will now repeat that procedure for a post-Newtonian (PN) 
effective potential which incorporates leading-order GR ef- 
fects for non-spinning, Schwarzschild black holes. Specifi- 
cally, we use the IPN harmonic coordinate Lagrangian pre- 
sented in |Blanchet| ( [20061 Eq- 174): 



^harm ^ Gmim2 _^ ITllvf 



2ri2 



C2 ( 



^'12 



+ 



mivf 



8 



+ 



Gm\m2 
ri2 



- I (^12 • 'yi)(m2 ■ 1^2) 



_L 3^2 



We define the effective potential as $ 



<=// 



_ 7 
4 

= K 



(42) 
{vi ■ 'U2))) 

harm 



where K represents the kinetic energy component of the La- 
grangian, i.e. those terms which depend only on velocities. 
This equation was derived for arbitrary mass-ratio systems, 
but here we identify the star as particle 1, the SMBH as 
particle 2, and have dropped all terms proportional to V2 or 
mi/m2. The Taylor expansion of $e/7 up to first PN order, 
around R — Rt, is given by 



GMbhR* 



1 + 



Rr 



2c2 Rt 



(43) 



where V is the star's orbital velocity at the tidal radius. 

From this equation, it is clear that GR corrections to 
the Newtonian potential will only matter for large, Mbh > 
lO^'M© SMBHs, with tidal radu close to or within the ISCO. 
However, all TDEs due to such black holes, or even more 
massive ones ( Kesden|2011 1 , will have debris energy spreads 
modified by GR around the ~ 2 level (although we caution 
that our PN approximation breaks down for tidal radii ap- 
proaching the ISCO). In this discussion we have neglected 
spin effects, but they may also play an important role for 
the subset of TDEs with Rp < -Risco- During completion of 
this paper, a more precise formalism was presented for esti- 



mating the GR corrections described in this section ( Kesden 
|2012) . The results indicate generally small (factors < 3) GR 
corrections which are maximized when the spread in energy 
of a spherical star is calculated near the ISCO, in qualita- 
tive agreement with our findings. We note however that |Kes"^ 
'den ('20I2I uses the older, less accurate approach to treating 
"frozen-in" debris energy, i.e. evaluating the spread in en- 
ergy at Rp rather than Rt. The primary difference between 
our approaches is that if debris energy freezes in at the tidal 
radius, GR corrections to Ae can only reach large, ~ 2 lev- 
els for Mbh > IO^'^Mq. TDEs around low-mass SMBHs wiU 
see negligible GR corrections to Ae, even if /3 is large. 



7 GRAVITATIONAL WAVES 

The parabolic motion of a star past orbital pericenter will 
produce low-frequency gravitational waves due to time vari- 
ation in the quadrupole moment of the star-SMBH system. 
For pericenters Rp < Rt, past work has indicated that such 
a signal could be marginally detectable with a LISA-like 
instrument ( Kobayashi et al. 2004 1 . Analogous work has 



focused on the inspiral of a white dwarf (WD) into an 
intermediate-mass black hole, where a similar signal could 



be generated by a violent disruption ( Rosswog et al.|[2009 
Haas et aL]|2012 1 , or a longer- lived G W signal could be ac 



companied by electromagnetic transients due to inspiral and 
stable mass transfer |Zalamea et al.|2010 ). An alternate, in- 
ternal source of GWs in TDEs comes from time variation 
of the star's own quadrupole moment during the Phase III 
vertical rebound, which was estimated in the past to gener- 
ate gravitational waves with strain h cc (Guillochon et 
al.|2009[ ) 



In this section, we present more detailed estimates of 
the "internal" GWs due to stellar pancaking and rebound, 
a process analogous to GW generation during core-collapse 
supernovae (CCSNe). As we shall demonstrate, TDE GWs 
are weakened relative to those in CCSNe due to lower col- 
lapse velocities and bounce accelerations, but are increased 
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due to the large degree of stellar asymmetry, and perhaps 
also by the correspondingly long lever arm of collapse in the 
quadrupole moment tensor. 

Specifically, we consider GW emission at the moment of 
maximum stellar compression, which we for now take to be 
synchronized throughout the star (but which will actually 
occur at different times for each point in the star, as seen 
in §5). The two polarization components of a GW signal, 
h+ and hx , can be read off of the transverse traceless GW 
strain 



'2G yTT 



(44) 



where d is the distance from the observer to the source. 



Jij — lij 



\5i-i5^^Iki is the reduced quadrupole moment 



tensor, jj^^ is a projection of J^j 
quadrupole moment tensor: 



and Iki is the standard 



d^fpx 



( 2x^ + 2xx 

2xy + xy + xij 
\^2xz + xz + xz 



(45) 



2xy + xy + xij 

2f + 2yy 
2yz + ijz + y'z 



2xz + ajz + x'z^ 
2yz + ijz + y'z 
2z^ + 2zz , 



To order of magnitude in the limit of fully synchronous ver- 
tical collapse, and neglecting the (weak) /3 dependence of all 
x and y terms, we then have 

dV /30 /30 



V/3 P 



where we have taken 7 = 5/3 (as we will for the remainder 
of this section) and approximated z ~ z/tc ~ /3^GMt/R^. 

If only diagonal terms are considered, the second deriva- 
tive of the quadrupole tensor will be oc /S'^. However, the 
extremely steep dependence of the off-diagonal Ixz and 
7yz terms indicates that viewing angles not closely aligned 
with z could in principle observe copious GW production. 
We note that physically, /xz ~ lO/yz because for the large /3 
where GW emission is relevant, y will be aligned with the 
star's longest principal axis. 

However, two degrees of symmetry present in this prob- 
lem will substantially reduce the magnitudes of Ixz and ly^. 
The free solutions indicate that to lowest order, tidally free- 
falling bodies should possess reflection symmetry about their 
in-plane principal axes. There is also an additional symme- 
try of reflection about the orbital plane. If these symmetries 
are exact, the off-diagonal terms in Eq.|45]will integrate to 0. 
The in-plane symmetries are broken when Rt/R ~ 1, i.e. for 
deeply plunging disruptions around low-mass SMBHs. The 
orbital plane reflection symmetry is more robust, and likely 
can only be broken by misalignment between the orbital 
plane and SMBH or stellar spin, which is beyond the scope 
of this paper. For the remainder of this section, we treat 
GW emission from off-diagonal terms as speculative, but 
the large magnitude of these terms in the integrand should 
motivate future work on disruptions of spinning stars, or 
TDES around spinning SMBHs. 

For /? < /?d, three-dimensional desynchronization is 
unimportant and the star collapses almost simultaneously, 
emitting GWs with a peak frequency of ~ 1/tc, which, using 
the calibration of Tc ~ 8.5/3^*r* from the afhne model and 



one-dimensional hydro simulations (§4), gives 



(46) 



where we have normalized m* = M^/Mq and r, = R^/Rq. 
Low mass stars have an easier time achieving high frequen- 
cies; if we use the relation oc M°'* for main sequence 
stars with M« ^ Mq, we find /gw ~ 10 Hz at /3 = 15, for 
M. = O.IMq. 

These frequencies are located on the far edge of the 
Advanced LIGO band, although with steep /3 dependence. 
Because the three-dimensional desynchronization discussed 
in §5 results in the leading edge of the star collapsing before 
the trailing edge (which lags by a time Atso), the GW sig- 
nal will be smeared out over a range of frequencies between 
l/AfsD and I/tc when /3 > /3d. This smears out the grav- 
itational wave emission by a factor AtsD/rc, giving us the 
strain estimates 



1 X 10"^^ Emlr~'^dul 




5/12 



(47) 



(48) 



Here dio is distance to the source normalized to 10 mega- 
parsecs, and for clarity (i.e. to separate diagonal and off- 
diagonal components of Jij into different polarization states) 
we have assumed a line of sight along the y axis so that 
dh+ — G(Jxx — Jzz)/c'^, and d/ix — 2GJxz/c'^- As in §5, 
/3d ~ 6 is the critical /3 value above which TDEs expe- 
rience significant three dimensional desynchronization. We 
have defined a parameter, H (which is most likely <C 1), to 
parametrize the unknown degree of reflection asymmetry in 
Phase HI of a TDE. 

The prospects for high frequency GW observation of 
TDEs involving main sequence stars appear dim, unless H > 
0.1. If we limit ourselves to GWs from J^^, then disruptions 
of abundant low mass stars have an easier time falling within 
the Advanced LIGO band, but produce too little strain to 
be detected; disruptions of solar-type stars produce a barely 
detectable strain at d ~ 1 Mpc, but will lie outside the 
Advanced LIGO band for /3 < 25. 

Tidal disruptions of white dwarfs by intermediate mass 
black holes appear more promising: a WD of mass 1 Mq 
and radius 6 x lO'^ m disrupted at /3 = 5 by a lO^M© IMBH 
will produce /i+ « 6 x 10~^* from a distance of 20 Mpc 
(neglecting all strain from off-diagonal terms in Jij). The 
peak emission frequency /ow ~ 60 Hz, but the existence 
of IMBHs is sufficiently uncertain that we do not attempt 
a rate estimate. We note that the level of emission seen at 
~ 50 Hz frequencies in full numerical relativity simulations 
of WD-IMBH disruptions (Haas et al. 20121 was approx- 
imately X 10'^ smaller than our prediction. This may not 
surprising, as the maximum density enhancement seen in 
these simulations is < 10, not ~ 200 as predicted by one- 
dimensional models (7 = 5/3, /3 « 10). It is likely that 
this discrepency is at least partially due to insufficient ver- 
tical resolution: after disruption, the smallest grid cell in 
these simulations is ~ 7?wd/40. Furthermore, Newtonian 
and pseudo-Newtonian simulations of /3 ~ 10 WD TDEs 
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( Rosswog et al. 2009 1 find maximum degrees of vertical com- 



pression ~ 100, in approximate agreement with the argu- 
ments presented in §4 and 5. 



8 OBSERVATIONAL IMPLICATIONS 

In this section, we discuss the observational implications of 
a revised Ae for optical transient searches. In general, we 
predict longer decay times but lower initial mass fallback 
rates than prior works ( Strubbe fc Quataert|2009| . 

Using the introduced index n, we rederive here a num- 
ber of important consequences of Eq. |4] in a parametrized 
way. The fallback time for the most tightly bound debris is 



itaii — 3.5 X 10 sec k 



3/2^-3n/2„,fl/2 



1 3/2 



(49) 



.0. 



where we have used the normalizations Ma 
Mbh/(1O'^M0), m. = M./Mq, and r. 
The rate of mass fallback evolves with time t as 

-5/3 

Mfaii ~ ( 7^ 1 , (50) 



t 



and is initially super-Eddington for disruptions of solar-type 
stars by SMBHs with Mbh < 10^'^Mq. Assuming a radia- 
tive efficiency < rj < 1, the peak (i.e. time of first pericen- 
ter return) mass fallback rate is given by 



A'/poak 



133 ^?_lfc=''''^/3='"'''M-^''^mJr-^''^ (51) 



implying that the maximum black hole mass that can un- 
dergo a phase of super-Eddington accretion is given by 

X = 2.6 X 10^ Mq rj^^klS^mt^^r-^. (52) 



Mb 



The mass fallback rate becomes sub-Eddington at a time 



tEdd = 6.6 X lO'^ sec ri!lk~'-^^^l3 



3/5 j,-3/5 ^-3n/5^,^-2/5^^1/5^3/5 



(53) 



although if t^dd < tfaii there is no super-Eddington accre- 
tion phase. Here we have set ij-i = rj/.l. Considering the 
dependences of Eqs. |49| |51| |52| and |53] on n we see that 
adopting Eq. [S] can have dramatic effects on TDEs from 
stars on deeply plunging (/3 > 3) orbits. In particular, us- 
ing the correct values of Ae reduces the peak mass fallback 
rate, decreasing the maximum SMBH mass that can produce 
a super-Eddington accretion phase. On the other hand, for 
TDEs with super-Eddington accretion, the duration of the 
super-Eddington phase is extended for (realistic) low-n val- 
ues. 

We plot the effect of n on the mass fallback rate in Fig. 
11 In previous literature (n = 2) a wide variety of mass 
fallback curves were possible, with high peaks and fast de- 
cay times accompanying large /3 values. If n = 0, however, 
the mass fallback rate is generally independent of /3. Un- 
der simplifying assumptions about the relationship of disk 
luminosity to M (often but inaccurately taken as L oc M; 
for complications see Lodato et al. (20091; Lodato & Rossi 



( 2011 1), the fallback timescale can be inferred by sufficiently 
long lightcurve observations. Alternatively, in the future it 
may be possible to measure tfaii directly, by measuring the 
delay between the onset of accretion and a prompt signal 
accompanying stellar disruption (either X-ray shock break- 
out or GWs). In either case, the j3 independence of tfaii will 



1000 




150 200 

t[d] 

Figure 11. Mass fallback curves (normalized by tlie Eddington 
fallback rate) versus time since disruption in days, for a variety 
of TDEs. SMBH masses of IO'^Mq, IO'^Mq, and 5 X lO'^ Mq are 
plotted as blue, orange, and purple curves, respectively. Likewise, 
[5 values of 1, 3, and 10 are plotted as solid, dashed, and dotted 
curves assuming n = 2, the power law index defined by Ae oc /3" - 
see Eq.|4] If n = 0, the solid curves represent all /3 values. Here we 
consider solar-type stars, and for simplicity set fc = 1 and r; = 0.1. 



simplify parameter extraction, in particular measurement of 

AfBH. 

Adopting n — will also alter the distribution of /3 in 
the TDEs detected by individual wide-field surveys, NdctiP). 
This is a quantity distinct from the distribution of the in- 
trinsic TDE rate, Ntde{I3), which scales as TVtde oc 
for two-body relaxation in the 



pinhole" regime (Brassart & 
Luminet|2008 1. Alternatively, if the dominant source of loss 
cone fueling is two-body relaxation in the "diffusion" regime, 
almost all TDEs will have P — 1; however, since most 
SMBHs are supplied with stars coming from the boundary 
between these regimes we will consider the "pinhole" regime 
for the remainder of this section (as it is the relaxational 
regime with nontrivial /3 dependence). 

If we first consider UV or soft X-ray surveys sensitive to 
the peak frequencies of TDE disk emission, then the peak 
luminosity L oc Afpcak oc implying a survey hori- 

zon rhor oc /J^"/" and a detection rate Ndet oc /g^^""*^/*. 
Optical detections of TDE disks will not be sensitive to 
the event's bolometric luminosity but rather to emission on 
the Rayleig h- Jeans tail, for which L oc M^''^ and Adct oc 
^(9n-i6)/i6 (jLodato & Rossi||201l|. However, both of these 



scaling relations for A'det assume a purely ffux-limited sur- 
vey; in the old picture of TDE energy spread (n = 2), high-/? 
events would be favored by their high ffux but disfavored by 
their shorter timescales of peak emission. If the timescale 
of peak emission, tfaii, is less than the survey cadence, tcad, 
then the probability of detection will be approximately re- 
duced by the factor ftaii/icad. This gives iVdot oc /3(3"-4)/4 
and iVdct oc _a(-i5"-i6)/i6 for X-ray and optical disk emis- 
sion, respectively. 

Although the details remain uncertain, several recent 



papers (Strubbe & Quataert 2009 Lodato & Rossi 20111 



have predicted that super-Eddington, radiation-driven out- 
ffows may dominate early emission from TDE accretion 
disks, particularly at long wavelengths. Using a simple black- 
body model with peak frequency fbb ( jLodato fc Rossi^201l1 , 
which predicts a peak luminosity L oc /S"'^"^* at frequencies 
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Scenario 




" = 2(^<1) 


n = 


Disk, 1/ l/bb 


7/2 


1/2 


-1 


Disk, u < i/bb 


1/8 


-23/8 


-1 


SE, f l/bb 


-2/5 


-17/5 


-7/10 


SE, U < l^bb 


81/16 


33/16 


-33/16 



Table 1. Here we display the scaling exponents s for the /3 de- 
pendence of the TDE detection rate, Ndetil^) ex /J". In the first 
column we describe the frequency and source of emission (SE in- 
dicates a super- Eddington outflow); in the second we give s for 
the standard scenario n = 2 and a high-cadence survey; in the 
third we again consider s in the standard scenario, but for a slow- 
cadence survey; in the final column we give s for our revised Ae, 
with n = 0, at any cadence. In this table t'bb refers to the "black- 
body" frequency of peak emission (the super-Eddington outflow 
is assumed to have a thermal spectrum; the disk emission is better 
modeled as a multicolor blackbody). 



< i^hb, or L oc /S^''^ for u > Uhb, we can repeat the above 
calculations. The P dependence of all different scenarios are 
presented in Table 1. With n — 0, all observational strate- 
gies are biased in favor of low /3 detections, including many 
strategies that once favored high (3 TDEs. 



9 DISCUSSION 

In this paper, we have analyzed the tidal disruption and free 
fall of a star in the context of "free solutions" to the Hills 
equations of the parabolic restricted three-body problem. 
The important conclusions of this work are the following: 

(i) During the tidal disruption of a star, debris energy 
freezes in at R = Rt, not R — Rp. 

(ii) Consequently, the spread in debris energy is smaller 
than in past analytic predictions. This will result in flares 
with longer fallback times, i.e. ones that decay more slowly 
but have smaller initial fallback rates. Fewer TDEs will drive 
powerful super-Eddington outflows than has been predicted 
in the past. 

(iii) The spread in debris energy is generally dominated 
by the freeze-in energy, although redistribution of the kinetic 
energy of vertical collapse to in-plane motions may result in 
slight variation in Ae for low /3(< 5). Rapidly spinning stars 
may see a stronger version of this effect at high j3 if their 
spins are misaligned with the orbital angular momentum 
vector. 

(iv) The leading order GR corrections to the frozen- 
out value of Ae are small, and generally negligible unless 
Rp ^ 67?g. For such TDEs (i.e. all TDEs for SMBH masses 
above IO^'^Mq) we have derived for the first time the IPN 
modifications to Ae. 

(v) The free solution model we have introduced is, in the 
limit of spherical and static initial conditions, an approxi- 
mate simplification of the affine model . However, the ability 
to include a range of nonspherical or dynamic initial condi- 
tions gives it a degree of flexibility not present in the affine 
model. Furthermore, the deformations to an initially spher- 
ical body in tidal free fall are not ellipsoids, as is assumed 
by the affine model. 

(vi) Gravitational waves are generated from variation in 
the internal quadrupole moment of a tidally disrupting star, 
reaching peak amplitude at the moment of maximum com- 



pression and bounce. GW emission is likely dominated by 
a single term in the quadrupole moment tensor, I^z oc P'^. 
For main sequence stars, these are unlikely to be detectable 
by ground-based GW interferometers. GWs from the tidal 
disruption of WDs by IMBHs are more promising targets, 
and for modest /3 values (~ 5) could be detectable to tens 
of megaparsecs. 

(vii) Gravitational wave emission from nondiagonal terms 
in the star's quadrupole moment could alter the previous 
conclusion, since Ixz,Iyz oc . However, in order for these 
terms to possess a nonzero pref actor, the refiection symme- 
try of the TDE about the orbital plane must be broken by 
either SMBH or stellar spin. Whether this can be done with- 
out desynchronizing the vertical collapse and weakening the 



The existing hydrodynamical literature did not until 
very recently support the first three of these conclusions. 
With a few exceptions, prior work has focused mainly on 
the common /3 = 1 events, for which Eqs. [2] and |3] are 
identical. In Laguna et al. (19931, the authors conducted 



SPH simulations in a static Schwarzschild background ge- 
ometry for /3 = {1, 5, 10}. They found that the fallback time 
scaled approximately as tfaii oc I3~^'^ and the velocity of 
unbound ejecta Va] oc Z?"'^, both of which imply the inter- 
mediate value of n = 1, for the power-law index defined in 
Eq. |4]as Ae oc However, the limited resolution of their 
simulations (7000 SPH particles) makes it unclear whether 
they would have possessed the midplane resolution to re- 
solve the phase of maximum compression; indeed, they find 
that pc oc a scaling well below analytic predictions 

as well as the higher-resolution one-dimensional simulations 
of BL08. 

It is also possible that our approximations have under- 
estimated the efficiency with which the bounce phase can 
redistribute the energy of vertical collapse to motions within 
the orbital plane, perhaps due to neglect of GR effects dur- 
ing phase III. Alternatively, the simulations of |Laguna et al.| 
( |1993[ ) may be altered by entropy production due to artificial 
viscosity, as suggested by the authors. 

More recent hydrodynamical simulations of high-/? dis- 
ruptions did not publish details on the spread of debris 
energy (Kobayeishi et al. 2004 Guillochon et al. 20091, al- 



though we note that Kobayashi et al. (20041, using a one- 



dimensional mesh code, found results in qualitative agree- 



ment with BL08. On the other hand, Guillochon et al. ( 2009 1 
found a vastly lower degree (pc ~ 4p*) of central compres- 
sion for a /3 = 7 event than what is predicted analytically, 
by the affine model, or by one-dimensional hydrodynamic 
simulations, with the difference attributed to a combination 
of insufficient midplane resolution and the physical, three- 
dimensional effects discussed in §5. However, we have shown 
that these multidimensional effects are unlikely to be large 
due to the proximity of the sonic point to the region of max- 
imum compression. 

During the completion of this paper, an independent nu- 



merical study ( Guillochon & Ramirez-Ruiz 2012 \ was posted 



which found an approximate constancy in Ae for 1 < /3 < 4, 
in the course of a thorough investigation of low-/3 TDEs. 
This result is broadly compatible with our paper, although 
further hydrodynamical simulations will be required to ex- 
plore high /3 values, and to test our prediction of desynchro- 
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nized collapse for rapidly rotating, misaligned stars (which 
could restore a /3 dependence to Ae, for deeply plunging dis- 



ruptions). Like Guillochon & Ramirez-Ruiz (2012 I, we argue 
that the assumption of frozen-in debris energy is made much 
more valid by assuming the energy freezes in when internal 
forces become negligible at R — Rt, not at later times. 

We note that the free solutions are of course an ap- 
proximation to the physical reality of tidal disruption, and 
neglect a number of important physical effects. Some areas 
in particular need of further clarification by hydrodynami- 
cal simulation are the velocity and shape perturbations in- 
duced on the star by the nonlinear hydrodynamics of tidal 
disruption at R ^ Rt. Although BL08 found these effects to 
generally be negligble for Phase II and III of a TDE, their 
simulations were one-dimensional and it is conceivable that 
in three dimensions the picture may change. The free so- 
lutions also fail to capture GR effects relevant for large /3 
values, although we have argued based on our own analy- 
sis and past literature that these effects do not qualitatively 
change our conclusions. Furthermore, we have not consid- 
ered the effects of SMBH spin, which if misaligned with the 
orbital plane may be able to desynchronize stellar collapse. 

Despite these limitations, the free solutions are a sim- 
ple yet powerful method for considering the tidal free fall 
of disrupted bodies. Although we have focused on the dis- 
ruption of main sequence stars by supermassive black holes, 
we stress that much of our analysis is general and applies 
equally well to other scenarios, such as tidal disruption of 
white dwarfs by intermediate mass black holes, or of planets 
or asteroids by more compact bodies. In all of these sce- 
narios, the decoupling of vertical from in-plane motion will 
strongly compress the disrupted object, often by many or- 
ders of magnitude. In future work, we aim to apply the free 
solutions to these alternative physical regimes, as well as 
to internal motions such as pre-disruption rotation or post- 
disruption shock formation. 
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APPENDIX A: PRINCIPAL AXES OF FREE 
SOLUTIONS 

Using the in-plane free solutions, we find that the initial 
phases of the principal axes at / = /c are given by 



cos 9c 



2 2^f + 4fe2 



where 



16 148 1312 

25^^ 25 ^ 25/3 5/?2 



352 

5/?2 

788 2208 



(Al) 



(A2) 



+ 



32 



k =\//3 



72 
25 



25 25/3 5/32 
936 1144 16 
25/3 25/32 " 5/33 
72 216 _ 176 \ 
25 25/3 ~ 5^^ 



are exact expressions. Note that this gives us 8 possible so- 
lutions for Sox- While 4 are spurious, the other 4 of these are 
valid, with each minimum (maximum) in ru possessing an 
equal minimum (maximum) 180° around the star's center of 
mass. In particular, for sin 60 > 0, 



COS&n 



and 



l + ^-^^'Tr7=f^ (A3) 



2^^ ^^'^2/fT4F' ^^"^^ 



Here the behavior of Smin and Smax is piecewise with 
respect to /3, due to zeros of j and k. We can describe this 
behavior for all /3 > 1 by setting 

r{l,2,2,l},/3< 1.073 
{p, q, P,Q}^ I {2, 1, 1, 2}, 1.073 s; /3 < 4.944 (A5) 
[{1,2,2,1}, 4.944 sC/3. 

To gain better intuition for the geometry of these prin- 
cipal axes, we can Taylor expand riong = »'H(Smax) and 
rshort = '■H(6'min) in the limit of large /3: 



Hong =g/3 
'"short —2/3 



1/2 



+ 22^-1/2^^(^-3/2) 

5 



-1/2 



23 , 



5/2n 



(A6) 
(A7) 



presented in Fig. [3] By expanding the numerator and de- 
nominator of tanVJc = — 2/H(6max)/a;H(6'max) we find 



tan-f/jc '■ 



16/3^/^ -I- 86^- 
8 + 185/3- 



-1/2 



(A8) 



Using trigonometric identities and the Keplerian expression 
tan To = v/?+ V/3 - 1 ~ 2/3^/^ (To is the angle between the 
X and a parabolic orbit's velocity vector at / = /c), we find 
the misalignment angle 

284 

32/33/2 + 180^1/2 + 185/3-1/2 • (A^) 

We have defined Vc as the (positive) angle between the center 
of mass velocity and the long principal axis at / = /o. It 
asymptotes to 0° as /3 goes to 00. 



These results are poorly convergent for /J < 10, but describe 
the exact solutions well above this threshold. We are now 
in a position to derive approximate formulae for the angles 
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